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Summary. — We give an overview of the current theory of collective modes in 
trapped atomic gases at finite temperatures, when the dynamics of the condensate 
and non-condensate must both be considered. A simple introduction is given to 
the quantum field formulation of the dynamics of an interacting Bose-condensed 
system, based on equations of motion for the condensate wavefunction and single- 
particle Green's functions for the non-condensate atoms. We discuss the nature 
of excitations in the mean-field collisionless region, including the Beliaev second- 
order approximation for the self-energies. We also sketch the derivation of coupled 
two-fluid hydrodynamic equations using a simple kinetic equation which includes 
collisions between condensate and non-condensate atoms. 
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1. Introduction 

There are several excellent review articles on trapped Bose-condensed atomic gases 
written at a relatively introductory level [1, 2], with emphasis on the dynamics of the 
condensate at T = 0. In the present lectures, I will concentrate on the interplay between 
the condensate and non-condensate components. This topic requires a more sophisti- 
cated analysis based on the concepts and methods of many body theory. My lectures 
will attempt to give a very basic introduction to this kind of approach. The primary 
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audience I have in mind are graduate students and postdocs coming from atomic and 
laser physics, rather than condensed matter theory. However, I hope the experts will also 
learn something from these lectures. In his lectures, Fetter [3] has given a detailed review 
of excitations in a trapped dilute Bose gas at T — 0. He shows that a convenient way of 
discussing these excitations is to start from the time-dependent Gross-Pitaevskii (GP) 
equation for the macroscopic wavefunction $(r, t). Linearizing around the equilibrium 
value (fr (r), one finds of the macroscopic wavefunction are given by the well-known Bo- 
goliubov coupled equations. For a uniform Bose-condensed gas, the excitation frequencies 
are those first discussed by Bogoliubov in 1947 [4]. 

At T = 0, one can assume that all the atoms are in the Bose condensate described by 
$(r, t). In contrast, in these lectures, my main topic will be to review our understanding 
of what excitations are in a Bose-condensed gas at finite temperatures, when there is 
a large number of atoms in the non-condensate (in trapped gases, the non-condensate 
is often referred to as the "thermal cloud"). I will make some contact with the ideas 
discussed in the lectures by Burnett [5] . 

The first half of these lectures (Sections 2-7) deals with excitations whose very ex- 
istence depends on self-consistent mean fields (of various kinds!), rather than on the col- 
lisions between atoms. In the standard language developed in condensed matter physics 
in the 1960's, this means the excitations are in the "collisionless region" . 

The second half of the lectures (Sections 8 - 12) deals with excitations in the collision- 
dominated hydrodynamic region. I review the two-fluid hydrodynamic equations such 
as given by Landau [6], generalized to include a trap. I give an explicit microscopic 
derivation of such two-fluid equations in a trapped Bose gas. This extends recent work 
of Zaremba, Griffin and Nikuni (ZGN) [7] to include the case when the condensate is not 
yet in local equilibrium with the non-condensate atoms. 

The type of questions I want to address in these lectures include: 

1. What is the difference between an elementary excitation and a collective mode? 

2. Can we isolate the dynamical role of the condensate on the nature of the excita- 
tions? 

3. At T — 0, with a pure condensate, excitations in a gas must be oscillations of 
the condensate. In contrast, above Tbec> the excitations are not related to a 
condensate. What happens to an excitation as we go from T = to T > Tbec- 
How does the excitation get rid of its "condensate" dressing? 

4. What is the essential physics behind the different mean-field theories of excitations 
which have been discussed in the recent literature: Gross-Pitaevskii, Hartree-Fock- 
Popov, Hartree-Fock-Bogoluibov? 

5. What is the physics behind the dreaded Beliaev second-order approximation? This 
is the first approximation which includes damping of the elementary excitations 
evenatT = [8]. 
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6. Why are the excitations and collective modes in a Bose-condensed system uniquely 
interesting, compared to all other many body systems? The key reason is, of course, 
that the condensate couples and hybridizes single-particle excitations with density 
fluctuations. Above Tbec? these two excitation branches are uncoupled. Once this 
is understood, one sees why more detailed and systematic experimental studies of 
excitations in trapped atomic gases at finite temperature are of great importance. 

To give a careful discussion of all these questions would need 10 lectures. In these 3 
lectures, I will give a speeded-up version. I will often use a uniform weakly interacting 
Bosc gas to illustrate the structure of the theory. While I will only sketch the math, I 
will still try to give a flavour of what is involved. The approach I will use to discuss 
these questions is based on the field-theoretic formulation of a Bose-condcnsed system of 
particles. As I review elsewhere in this volume [9], this powerful formalism was introduced 
by Beliaev [10] in 1957 and extensively developed in the Golden Period: 1958-1965. I will 
introduce this formalism in a very schematic manner - but even experimentalists will find 
it useful to know some of the "language" used in this approach. There are other methods 
to deal with collective modes in Bose gases but they are not as useful at isolating the 
dynamical role of the condensate, or dealing with finite temperatures. 

While I will always have trapped atomic Bose gases in mind, much of the general 
theory [11, 12] is valid for any Bose-condensed fluid (gas or liquid). Thus I will often 
make references to superfluid He, pointing out similarities with Bose gases. 

2. Elementary excitations and density fluctuations in normal systems 

We work with quantum field operators : 

V>(r) = destroys an atom at r 

(1) V> + (r) = creates an atom at r. 

These operators satisfy the usual Bose commutation relations, such as [f/>(r), -0+(r')] = 
<5(r — r'). Of course, if we have several hyperfine atomic states, then we have different 
field operators V> a (r), where a is the hyperfine state label (the analogue of a spin label). 
All observables can be written in terms of these quantum field operators. Two important 
examples are the local density operator: 

(2) n(r) = Vi + (rMr) 
and the Hamiltonian 



(3) 



H = J dr^+(r) -X! + [/ cxt (r)- M $(r) + J dr^+(r)^+(r)^(r)^(r). 



In these lectures, the two-particle interaction will always be approximated by a s-wave 
scattering length, appropriate to a dilute Bose gas at very low temperatures [5, 13]. Thus 
v(r — r') = gS(r — r'), with g = (we set K = 1 throughout this article). 
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What we want to calculate are various kinds of correlation functions involving different 
local operators: 

(h(r,t)h(r',t')) ~Xnn(rt,rV) 

= Xrm(l, 1') = density response function 
(j a {v,t)h{r't')) -X^KrY) 

(4) = Xj a jpO-i 1') = current response function. 

In these tutorial lectures, I will not bother to distinguish between the various kinds of 
correlation functions (time-ordered, retarded, etc) or even the difference between correla- 
tion functions and response functions. These are discussed in all the standard textbooks 
on many body theory [12, 14] but are only important when one is doing systematic cal- 
culations. A convenient summary of the formalism is given in a recent review article on 
homogeneous weakly interacting Bose gases [15]. 
The density response function 

(5) X n„(l,l') = ^ + (1)^(1)V' + (1'W1')> 

involves four quantum field operators and is an example of a two-particle Green's function 
6*2(1, 1')- One can measure x„ n (l,l') by coupling a weak external field to the local 
density in the system 

(6) H'(t) = J drV ex (v,t)n(r) = J dlV cx (l)n(l) . 

Linear response theory gives 

5n(l) = (n(r)>, - (n(r)) e , 

(7) = J dl' Xnn (l,l')V ex (l') + ... 

for the density response (see Chapters 5 and 9 of Rcf.[12] and Chapter 2 of Ref.[14]). For 
a uniform system with V ex (r, t) = V q ^ e*( q ' r_wt ) , we have 

Xnn(l,l') =Xnn(l-l') 

= Xnn(r-r',t-t'). 

Fourier transforming (r — r') — > q and (t — t') — > ui, the linear response equation in (7) 
reduces to 



(8) 



Sn(q,u) = Xnn(<l,v)V qiU1 . 



THEORY OF EXCITATIONS AT FINITE TEMPERATURES 



5 



If 

Xnni^L,^) ~ u-e nas a P°^ e &t w — E q , (for further details, see Section 3) then when 
u and q of the external potential satisfy cj — E q , Sn(q,uj) as given by (8) can be very 
large even though V q ^ is small. We note that 

(9) Xnn(q, W) l — -> Xrm(q, t) ~ e~ iE ^ 

uj — £j q 

and thus clearly the pole at lo = E q is the signature of an oscillating density fluctuation. 
The basic correlation function in the field-theoretic approach is given by 

(10) TO+a'D-dlM'). 

This single-particle Green's function G\ involves two quantum field operators. It de- 
scribes creating an atom at 1' = r',t', let it propagate through the system to 1 = r, t and 
then destroying the atom. All other higher-order correlation functions such as C?2(l, 1') 
can be constructed out of combinations of Gi(l, 1'). For example, the lowest order con- 
tributions to the density response function in (5) are 

x„„(i,i') = ^ + (iWi^ + (i'^(i')> 

~ $+(l$(l))$+(l'$(l')> (n(l))<n(l')> 

+ (^(1)^(1'))(^(1>(1)) - d(l, l')G!(l', 1) 

(11) +(Vi+(l^+(l'))(Vi(l)Vi(l')) + ... 

As we discuss in Section 4, the terms in the last line of (11) vanish in a normal Bose 
system but are finite for T < Tbec- The poles of G\ correspond to what are called 
elementary excitations (or quasiparticlcs) . In a uniform system, we have 

Gi(l,l')->Gi(q,o/) 



a 2 

(12) Gi(q, t) ~ e- lE « l ; in a free gas, we have E s q p = ^— = e g . 

One can show [12] that these single-particle excitations determine the thermodynamic 
properties of interacting systems. However, it is very difficult to directly measure the 
spectrum of Gi(l, 1') since one needs an external field which couples to $(1), i e 7 an atom 
reservoir. Later we will see that what makes a Bose-condensed system unique is that we 
can easily access Gi(l, 1') as a result of the effects of the Bose condensate. 

Finally we introduce the key idea of a single-particle self-energy through Dyson's 
equation: 



(13) 



Gi — Go + GqSGi, 
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where G\ is the interacting single-particle Green's function, Go is the non-interacting 
single-particle Green's function and all effects of the two-particle interactions are con- 
tained in the self-energy function E. In a uniform system, we can Fourier transform this 
Dyson equation to give 

(14) G (q,oj) = G {q,uj) + G (g,w)I!(g,w)Gi(gr,w), 
where Go{q, oj) = . This is easily solved to give 

(15) Gi = p-^— - y 

Thus we see that G\ may have a single-particle pole at the quasiparticle energy 

(16) E s / = e q + X(q,E°v). 

In general, u>) — E# + iEj, where Ej describes the damping of the single-particle 
excitations. 

Field-theoretic calculations [11, 12] involve a systemmatic (diagrammatic) procedure 
to calculate Ti(q,u>) and from this to obtain G\(q,u)). We note that the self-energy E is 
by definition highly non-perturbative. As an illustration, let us consider the self-energy 
to first order in g. This Hartree-Fock approximation is shown in Fig.l. In our s-wave 
approximation, the total self-energy is simply 



(17) E^f = gn + gn = 2gn, 
and therefore (15) gives 

(18) G(q,w) = ; ] 



[e q + 2gn) ' 

Hence the normal HF excitation energy has the dispersion relation 

(19) Er = £+2gn. 

3. Density fluctuation spectrum in the mean-field approximation 

The simplest approximation for the density response function introduced in Section 2 
is [11, 12] 

(20) Xnn(h 1') => Xnnih l') = G^l, l')Gl(l', 1). 
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In a uniform Bose system, the Fourier transform of this gives 

(21) £ n (q,u>) ~ fdk fdLJ! [ dcj 2 ^ N ° {UJl) - N ° {U,2) ^ A(k, ^)A(k - q, u 2 ). 

J J J Ul — u 2 — u 

Here N°(u) = (e' 3 " — 1) _1 is the Bose distribution function and A(k,u>) = single-particle 
spectral density ~ 2ImGi(k, u> + i0 + ). If we use AnF(k, w) ~ 2ir8(oj — Ej^ F ) as given in 
(18), we find (21) reduces to 

(22) Xnn («, W) ~ J dk _ E H_F ) _ LJ ■ 

We note that x«„ nas a continuum of poles given by w = (-E^ F — E^). It is easy 
to understand the physics which gives rise to this "ideal-gas" spectrum. The Fourier 
transform of the local density operator in (2) is given by (for a uniform system) 

(23) h q = ^2a+a k - q 

k 

Clearly h q creates a "particle-hole" density fluctuation with the following features: 

change in energy : E^ — E^_ q = lo 

(24) change in momentum : k — (k — q) = q. 



The spectrum (22) describes a broad incoherent superposition of particle-hole states [14] 
and is not a true collective mode, such as we discuss next. 

The mean field approximation (MFA) for the density response is also called other 
names: SCF (self-consistent field), RPA (randon phase approximation), but all involve 
the same physics. The MFA was introduced by Bohm and Pines in the period 1951-1953, 
work which had a pivotal effect in our understanding collective effects in all many-particle 
systems. We recall the linear response expression in (7), where x„„ is the full density 
response function for interacting Bose gas. If we introduce the self-consistent Hartree 
mean-field: 

(25) -^HartreeW = / dl ' v ( l - I'W) = g5n(l), 
then we can approximate the linear response equation in (7) by 

(26) Sn(l) = J dl' X ° nn (l, l')[Kx(l') + ^Hartree( 1 ')]- 

The system is assumed to respond as if the atoms propagate independently (as described 
by x ) but are moving in an effective field V e ff (!'). In a uniform system, the Fourier 
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transform of (26) is 

(27) Sn(q,cj) = Xnn(q,u)[V qu + g6n(q,uj)}, 
which gives the well-known MFA expression for Xnn [12, 14]: 

(28) Xnn{q,u) — q p r- 

This result for Xnn(q, w ) ma Y have new poles given by zeros of the denominator, 

(29) l-g X lM^)=0- 

Assuming that it is distinguisablc from the incoherent ideal gas density fluctuation spec- 
trum given by (22), this pole at w — E^ oU is called a zero sound mode (a plasmon in 
charged systems). This language was first introduced in 1957 by Landau in Fermi sys- 
tems [14] but the concept is generally applicable in any interacting many body system. 
Physically, it is clear that zero sound is a "collisionless" density oscillation arising from 
dynamic self-consistent mean fields. 

We make a few comments about such "zero sound" collective modes: 

1. E%° u and E s v (the poles of Gi(q, u)) are both states of an interacting many-particle 
system. However the collective mode disappears if there are no interactions, while 
the single-particle excitations still exist in a non-interacting gas. 

2. Because of the low density, dynamic mean fields in normal systems are too weak 
to allow the existence of a well-defined (weakly-damped) zero sound mode in the 
Bose gases of current interest (for T > Tbec)- If gx«n(?! w ) ^ 1; we can then 
approximate Xnn(q,u) ~ x°„(<Z,^)- The situation is quite different for T < Tbec, 
when a coherent mean- field due to the condensate is present. 

3. One also expects a collective pole to appear in Xrm(q, w) in the collision-dominated 
hydrodynamic region. However this sound wave pole is induced by rapid collisions 
producing local equilibrium, as we discuss in Sections 8 - 12. Ordinary sound is 
not described by the MFA discussed above, ie, it is not the result of dynamic mean 
fields. 



4. Green's function formulation of excitations in a Bose-condensed system 

All the formalism we have been discussing for T > Tbec can be extended in a natural 
way to Bose-condensed systems [10, 11, 12] making use of the fundamental decomposition 
which separates out the condensate and non-condensate parts of the quantum fields: 



(30) 



^(r) = <^(r))+^(r). 
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Here the average is over a restricted ensemble [16, 17] consistent with (ip) ^ 0. The most 
profound way of doing this is to add a symmetry-breaking field 



(31) H sb (t) = / drfo(r,t)^+(r)+H.C.] 



and work with Htot — H 'system + H s b, taking the limit r\ —* at the end. This gives the 
system a "hunting license" [16] to have finite value of (tp) and one finds 



(32) 



for T > Tbec ; {i>)sb ^ for T < T BEC . 



A key point is that if {tp} s b ^ 0, then a direct consequence is the existence of "anomalous' 
or "off-diagonal" propagators [10] 



(33) 



U+W+(i')Wo. 



These describe new condensate-induced correlations between non-condensate atoms at 
different space-time points 1 and V . In a sense, these anomalous correlation functions 
are as important as the macroscopic wavefunction $(1). We also note that m(l) = 
(ip(l)ip{l)) s b is the "pair" function that Burnett [5] discusses in detail, using a different 
formalism. ^From now on, we leave the symmetry-breaking label on the averages implicit. 

In addition, if (ip) ^ 0, we find that correlation functions involving three non- 
condensate field operators can be finite 



In particular, (n(l)^(l')) ^ 0. This describes the condensate-induced coupling of non- 
condensate density fluctuations n — ip + ip and the single-particle field fluctuations (see 
Section 7). 

Clearly one has to work with a 2 x 2 matrix single-particle propagator G\ when 
(t/j) 7^ 0, namely 



where we have introduced a 2 x 2 matrix Green's function for the non-condensate atoms 



(34) 



Mlty(2M3)) ± o. 



(35) 




$(i)iHi')> (VKiMi')) \ 

$+(i$+(i')) $ + (i$(i')> J 

$(1)$*(1') + Gn $(1)$(1') + Gi2 
$*(1)**(1') +G 2 i $*(!)$(!') + G 22 



) 



(36) 
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Beliaev [10] effectively showed (in modern matrix notation) 

(37) G a fi = GoSa/3 + G Y, a sG$f3, 

where we use the standard convention that repeated indices (S — 1,2) are summed over. 
This is the famous Dyson-Beliaev equation, involving a 2 x 2 matrix self-energy T, a s- 
Clearly all components G a p will share the same single-particle excitation spectrum. In 
a uniform system, we have 

(38) G a0 (q,u)) = G (q,Ld)S a p + G o (q,ui)T, aS (q,ui)G S 0(q,u)), 

which is a set of linear algebraic equations which are easy to solve for Gn(~ G22) and 
Gi2(~ G21). 

We will now use this Beliaev formalism to discuss various simple approximations for 
the self-energies S Q/ a. The interaction energy in (3), namely 

(39) V int = \g I rfr^+(r)^+(r)^(r)^(r), 



splits into various distinct contributions when we use (30) to separate out the condensate 
parts (see Fig. 2). At T = 0, we can ignore the V3 and V4 contributions because so few 
atoms are in non-condensate. This is the famous Bogoliubov approximation [4], and it is 
equivalent to the linearized GP theory. In Beliaev language, the Bogoliubov self-energies 
are shown in Fig. 3. For a uniform gas, these give [10, 12] 

ri , \ ui + € p + n c g 

Gu(p,w) - 



- [ep + 2e p n c g] 

(40) Gi 2 (p,w) 



^ 2 _ [ e l + 2e p n c9\ 

These Bogoliubov single-particle Green's functions contain the same physics as discussed 
in Fetter's lectures [3]. They clearly have poles at the frequencies u> = ±E P , where 

(41) E p = (e 2 p + 2e p n c g) 1 / 2 . 

At low p, the single-particle excitation is phonon-like E p = csogP, with the Bogoliubov 
phonon velocity 

(«) 4. = 

It turns out that this simple T = Bogoliubov theory already exhibits most of the 
structure which will always arise in Bose-condcnscd systems. This is why it plays the 
role of the "H-Atom" in discussions of Bose-condensed fluids [11]. 
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In this approximation, as we have stressed, the starting point assumes that all atoms 
are in the condensate and hence n c = n. However, we can still use Gu(p,w) in (40) to 
estimate the depletion. This makes use of the relation 

h = n — n c ~ Gn(l, 1) 

(43) = J dqj dwN°(Lu- ^)i n (q,u;), 

which is an exact formula for n (see Ch.4 of Ref. [11]) in terms of the so-called single- 
particle spectral density [see (21)]. The spectral density given by the Bogoliubov approx- 
imation above is 

(44) A Bog {q, w) = u\5{u) - E g ) - v 2 5(u> + E q ), 

showing the characteristic negative energy pole at u> = —E q . This last feature emphasizes 
that creating an atom in a Bosc-condensed system involves a coherent weighted combi- 
nation of creating an excitation and destroying an excitation. Calculating the depletion, 
one finds 



n 



dq J dujN°{uj)[u 2 q 5{uj - E q ) - v 2 S(lo + E q ) 
= J dq{N\E q )u 2 q -N\~E q )vl} 
dd{v 2 q +N°(E q )(u 2 q + v 2 q )} 

3\ 1/2 



(45) = -n u 

a result first obtained by Bogoliubov [4] and reviewed in detail by Fetter [3]. We have 
made use of the key identity which the Bose distribution satisfies, 

(46) N°{-E) = -[N°(E) + l}. 

We see from (45) that the T = depletion of the condensate formally arises from the 
negative energy pole in Gn(q,w). 

It is easy to evaluate the density response function Y° n (g,o>) in (20) using the Bo- 
goliubov spectral density in (44) . The negative energy poles are seen to give rise to new 
poles at uj — Ek + Ek- q which have finite weight even at T = 0. These are in addition 
to the "normal" particle-hole poles of \nn at u — Ek — E/-- q , as given by (22). 



5. Beyond the Bogoliubov approximation: classification in terms of self- 
energies 

At finite T, a lot of atoms are thermally excited out of the condensate in a dilute Bose 
gas. These "normal" self-energies must be included. In the Bcliaev self-energy formalism, 
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this was first done by Popov in 1965 [18]. We have to add the ordinary Hartree-Fock self- 
energies, as shown in Fig. 4. This gives the first-order Popov approximation (the excited 
atoms are treated as as ideal Bose gas, and thus the approximation is not self-consistent) 

Y,^° pov (p, w) = 2n c g + 2n a g 
[ l) £f 2 opo >,^) =n c g. 

These self-energies lead to [compare with (40)] 
(48) G n (p, e>) = ^ 2 P _ E2 ; G12 (p, w) = ^ 2 _ g2 , 



with 

(49) A = /ijjp — 2n° 3 ; £ p 2 + 2e p A, 
where 

(50) mP = S n (p = 0, u = 0) - Ei 2 (p = 0, w = 0). 

The "chemical potential" /inp as defined in (50) here was introduced by Hugcnholtz 
and Pines [19] at T = (and generalized to T ^ later [17]). HP showed that if 
G a p(p,u}) had a gapless excitation spectrum in the p — > 0,w — > limit, then the true 
chemical potential /x must satisfy fi — finp- In the Popov approximation, calculation 
gives [20] 

(51) Mhp = 2n c g + 2h°g - n c g = n c g + 2h°g, 

and thus A = n c g, where n c (T) is the ideal gas condensate density at temperature T. 
The generalized GP equation determines the chemical potential, which we find convenient 
to denote by \x c . In the Popov approximation, this equation is given by 

(52) j_Z_ + C / ext+5n o + 2fl n j$ = /ic$, 

and hence in a uniform gas, we see that /j, c = /Uhp- In the Popov approximation, the 
single-particle Green's functions G a ,p in (48) have poles given by 

(53) E p = [e 2 p + 2e p n c (T)g] 1 / 2 . 

This is formally the same as the T = Bogoliubov excitation frequency, except that now 
n c (T) is temperature-dependent. For a more detailed discussion, see Chapter 3 of the 
recent review by Shi and Griffin [15]. 

The full Hartree-Fock Bogoliubov (HFB) approximation involves calculating Y* aj 3 self- 
consistently using the complete matrix 2x2 Beliaev propagator, as shown in Fig. 5. 
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We note that above Tbec, the HFB reduces to the usual self-consistent Hartree-Fock 
approximation. Moreover, the HFB is the best single-particle approximation for a Bose- 
condensed system, in a variational sense. Within this limitation, the HFB will give the 
best single-particle approximation for the thermodynamic properties. 

However, there is one "bad" aspect of the full HFB which was noticed in the 1960's 
and which has been the subject of a several recent theoretical papers on Bose-condensed 
gases. Namely, the HFB docs not obey the Hugenholtz-Pines (HP) theorem. In a uniform 
gas, this means that the single-particle excitations will have an energy gap in the long 
wavelength limit q — > 0. It is easy to check this by calculating the two chemical potentials 
introduced earlier: 

,uhp = 2(n c + n)g - (n c + rh)g 
(54) = g(n + h — rh), using n = n c + n. 

The static equation of motion for <E> = ^fn~ c using the HFB approximation gives 

\fn c = [i c \fn c — ► n c = g(n + h + rh). 

One "solution" of this problem is to leave out rh (which is the source of the problem) 
but keep n. This corresponds to the self-consistent HFP(opov) approximation [20, 21, 
1] which, while approximate, has several nice features. It gives the correct excitation 
spectrum for both T — > and T > Tbec- Moreover, the spectrum is gapless at all 
temperatures. 

The origin of this "problem" with the HFB is clear. The HFB keeps all self-energies 
which are first order in the interaction g. However, by computing these self-energy 
diagrams using self-consistent G Q( a propagators, one is clearly bringing in terms to all 
orders in g. This brings one into dangerous territory! Moreover, one easily can check 
that rh must be at least of order g. Thus the gin contribution to £12 is at least of 0(g 2 ). 
This suggests that to "fix-up" the HFB as a theory of excitations, we have to include all 
self-energy contributions to at least second order in g. 

Moreover, it is not obvious that it is consistent to calculate n self-consistently [21] and 
ignore rh, since the lowest order interaction contributions to both quantities are found to 
be of the same order. Indeed, one can show for a uniform gas that keeping the lowest order 
asymptolic correction to the ideal gas result gives h = n cr — rh, where n cr = 7i(T/Tbec) 5 
is the critical density of an ideal gas at temperature T [22, 15]. However, we note that 
an expansion around the ideal gas results is not necessarily a good approximation in 
trapped gases (see also the discussion of the results in (106). At T = 0, the Popov 
approximation is a poor guide for corrections to the Bogoliubov approximation since a 
calculation analogous to (45) gives rh — 3h [15]. 

So, finally, we come to the work of Beliaev [8] who, in 1957, evaluated the E Q( a self- 
energies keeping all contributions up to order g 2 . This was originally done at T = 0. 
For finite temperatures T <~ Tbec> Bcliaev's work has been recently extended by several 



(55) 



2m 



+ £4xt + gn c + 2gh + grh 
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authors [18, 15, 22]. When one includes the second-order diagrams for Y, a p, we note that 
one must be careful to also treat the HFB first-order diagrams correctly to order g 2 since 
there are many terms which cancel. 

As expected, the Beliaev-type second order calculation cures all the problems of the 
first-order HFB approximation. We make a few comments on the results for a uniform 
Bose gas. The self-energies are found to be 



(56) 



En(g,w) = 2n c g R + 2n^g + A(q,uj)g 2 
£12(9, w) = n c g R + m (1) g + B(q,u)g 2 , 



with the renormalized interaction defined by 



(57) 
and 
(58) 



9r = g 



1 + 9 



r dv 
J W 



dk 1 



rn 



(i) 



-n c g 



dk 



) 3 2e k 



2N°(E k ) 



2E k 



Here the anomalous and normal densities, rh^ = Gi2(l, 1) and = Gn(l,l), are 
computed using the HFP approximation for G a /3- gn is the second-order approximation 
to the ^-matrix, as shown in Fig. 6. Using these results to calculate the chemical potential, 
we find (for q, lo — > 0) 



(59) 



/i H p = 2n c g R + 2h (1) g - n c g R - m (1) .g + (A - B)g 2 
= n c g R + 2h^g + m^g. 



This last result gives precisely the HFB result for y(i c , and hence the HP relation is 
satisfied. In the last line of (59), we have used the key result that for q, uj — > 0, (A—B)g 2 = 
2m( 1 )<7, as proven at all temperatures by Talbot and Griffin [23]. 

We note that one can rewrite the expression in (59) in the alternative form 



(60) 



M = 9 



n c + + 



where we have introduced a "renormalized" anomalous density 



(61) 



-n c g 



n c g 

dk 

(2^)3 



f dk 1 

J j2^Ye 



2e k 
2N°(E k ) 
2E k 



1 

2Tk 
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In contrast to in (58), rh^ has no ultraviolet divergence from large k contributions. 
Proukakis, Burnett and coworkers [5] have pointed out that the self-consistent ladder 
diagram approximation for the t— matrix (see Ref.[24] and Ch.4 of Ref.[15]) can be ex- 
pressed in terms of m/j. In a uniform gas, this relation is simply t = g(l + rhji/n c ). 
This suggests that the effect of m can be partially included by working in terms of the 
self-consistent t— matrix [5]. 

6. Physics of the Beliaev approximation 

We have seen in Section 5 that Beliaev's second-order calculation satisfies the HP 
theorem. It has another feature whose significance was only realized in the early 1960's. 
Using the above T = Beliaev expression (60) for the chemical potential fj,, one finds at 
T = (a result first found by Lee and Yang using a different approach) 



(62) 



MBci = ng 



3 \ 7T 



1/2- 



where we have used [see (45)] 
(63) 



n = n c + n c - 



8 / n r a c 



n c + n. 



3 V it 

One also finds that the phonon pole of Gf j' has a velocity given by 

„3\ 1/21 



(64) 



c Bcl 



4-Kna 



1 + 16 



again expressed in terms of the total density n. One may easily check that 



(65) 



2 n dfi(n) 



that is, the Beliaev phonon excitation at T — has a velocity precisely equal to the 
compressional sound velocity. 

This last result raises the basic question: why does the elementary excitation (pole of 
G a p) have a velocity corresponding to that of a density fluctuation (pole of Xrm)? Later 
this strange identity was proven to be correct at T — to all orders of perturbation 
theory by Gavoret and Nozieres [25]. To be precise, at T — 0, in any Bose-condcnscd 
fluid (liquid or gas), GN showed (for q,ui — > 0) 



(66) 



u A — c z q z 



uj 2 — c 2 q 2 
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where c 2 is given by the compressibility as in (65). We discuss this equivalence in Sec- 
tion 7. 

Let us consider the full time-dependent HFB equation of motion for $(r, t) = (^>(r)) t , 
given by 



(67) i 



d®(r,t) 
dt 



2 



•— + U ext (r) + gn c (r, t) + 2gn(r, t) + gm(r,t) 



$(r,t). 



There are various approximations to this key equation which have been used in recent 
discussions of trapped Bose gases: 

• GP: ignore h and m completely. This is valid at T = [1, 26]. 

• Static HFP: Keep n c (r,t) but set n(r,t) = n (r) and m(r,t) = [21, 27]. 

• Static HFB: Keep n c (r, t) but set n(r, t) — h (r) and m(r,t) — m (r) [17, 20]. 
As we have noted, this produces an energy gap in the excitation spectrum ( since 
mp ^ Mc)- 

• Dynamic HFB: Treat all dynamic mean-fields due to n c (r, t), h(r,t) and m(r, t) 
on an equal basis in a generalized mean-field calculation of the density response 
function x„ n (l,l'). One finds [28] that the poles of Xnn given by this kind of 
calculation (xnn ~ U i-E i ) are identical to the poles of G a p given by the second- 
order Beliaev approximation (G a p ~ ^-e'* )• 

Recent work on excitations at T ^=0 in trapped gases is easily understood in terms 
of such a generalized linear response calculation of Xnn- For example, Minguzzi and 
Tosi [29] only keep fluctuations in Sn c and Sn but dm is completely ignored (dynamic 
HFPopov); in contrast, Giorgini [30] keeps fluctuations in both 5n and 8m but only those 
induced by the condensate fluctuations 8n c . 

Let us briefly discuss the excitation frequencies in trapped gases at finite tempera- 
tures [5, 31]. For T < 0.4 T BE c, GP theory (T = 0) predictions [26] are in excellent 
agreement with the JILA [32] and MIT [33] experiments. For T > 0.6 T B ec, however, 
there are significant differences between most current theoretical results and the mea- 
sured excitation energies [21, 27]. The problem clearly is that mean- field theories must 
include the collective dynamics of both the condensate and the non-condensate. The 
Beliaev treatment does this, as do the papers [29, 30] mentioned above, as well as work 
based on a collisionlcss kinetic equation [7, 34] . These improved treatments can correctly 
obtain the in-phase rigid oscillation of the condensate and non-condensate equilibrium 
profiles (the so-called Kohn mode), as discussed in Section 11. 

7. Why are excitations in a Bose fluid so interesting? 



We have seen in Section 6 that at T = 0, the elementary excitations have the same 
spectrum as density fluctuations. In fact, one can prove this equivalence to be a direct 
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consequence of the Bose broken symmetry, $(r) = (ip) ^ 0. Thus, this equivalence is valid 
at all T < Tbec- Of course, as soon as T ^ 0, the phonon velocity is no longer simply 
given by the compressional sound velocity (65) as found at T — 0. This fundamental 
equivalence was only proven and understood in the early 1960's (Gavoret and Nozieres 
[25], Hohenberg and Martin [17], Bogoliubov [16]). It was implicitly assumed by Landau 
[6] in his 1941 paper and later by Feynman [35] in his 1953-54 papers on excitations. Of 
course, the clear distinction between elementary excitations and density fluctuations in 
many-body systems was only clarified in the period starting from 1957. Only then were 
Bose fluids realized to be very special [16, 17, 25]. Today, we understand this key hidden 
assumption of the work of Landau and of Feynman has its microscopic basis in $(r) ^ 0. 

Let us try to understand the origin of this equivalence of the field fluctuation and 
density fluctuation spectrum. The density operator in a Bosc-condcnsed system can be 
decomposed using (30) as follows: 

n(r) = i/; + (r)^(r) 

(68) = |$ (r)| 2 + $ (r)^ + (r) + &1(t)$(t) + ^ + (r)^(r), 

where the first three terms correspond to n c (r) = n c0 (r) + Sn c (r). Clearly density 
fluctuations in a Bose-condensed system (ie, $o 0) have a direct coupling to the single- 
particle field fluctuations, due to the possibility of atoms coming in or out of condensate 
reservoir. In momentum space, we recall that (68) is equivalent to 

(69) h q = ^2 o-t^k-q = \aa\ 2 + (a^ a— q + a a+) + h q . 

k 

Making the same decomposition as in (68) for the density response function (5) gives 

Xnn(l.l') = (n(l)n(l')) 

= (h(l))(h(l')) + (5n c (l)5n c (l')) 

+ (6n c (l)5n(l')) + (5n(l)5n c (l')) 

(70) + (<5n(l)<5n(l')>. 

We note that in the Bogoliubov approximation, we have Sn — and the equivalence of 
Xnn and G a @ spectrum is then trivial (see also [3]) 

(71) Xnn(q,u) Xnn(q,u) = 7l c0 ^ G^q, u). 

a.,0 

To illustrate how the hybridization of single-particle and density fluctuation arises, 
we discuss the time-dependent HFP approximation on which Ref. [29] is based. For 
simplicity, we consider the uniform gas case. In the MFA of Section 3, we have 



(72) 



Sh = x [V cxt + 2g{Sn c + Sn)}, 
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and hence for V e 



ext 







(73) 



Sn = 



X°2g5n c 



Using this result in the HFP equation of motion (67), ie, with 8fh — 0, we find 



(74) 



(2g5h + g5n c ) = g 



( 



1 + 2.9X' 
1 - 2<7X ( 



) 



Sn c = g'Sn c . 



In this approximation, which only keeps the Sh fluctuations induced by the condensate, 
the equation of motion for <5$(r, t) reduces to the GP equation, but with g — > g' '. Thus 
the characteristic poles will be determined by the zeroes of 



This is the same denominator (up to a common factor) exhibited by the response func- 
tions obtained by Minguzzi and Tosi [29] . This calculation illustrates how the Bogoliubov 
single-particle mode and the zero sound density fluctuation are coupled and hybridized. 
More generally, one can prove that (at all T): 



where A a (l, 1') is a Bose broken-symmetry vertex function which vanishes if n c o = 0. 
Using the "dielectric formalism" developed in the early 1970's (for a review, see Ch.5 of 
Ref. [11]), one can prove that the self-energies S Q /3 are such that \nn and G a p have the 
same poles — and that there are none specifically associated with \nn • The advantage of 
this diagrammatic formalism [36, 37, 38, 39, 40] 

is that one manifestly sees that the poles of G a p and Xnn are identical, within a given 
approximation for the "building blocks" of the dielectric formalism. We now see why we 
can study elementary excitations from the density fluctuation spectrum when $o ^ 0. 
To directly measure Gi(l, 1'), we need an external perturbation which can change the 
number of particles in the system. BEC indirectly provides us with such a probe! 

This equivalence of the single-particle excitations with density fluctuations lies at the 
heart of the phenomenon of superfluidity It essentially restricts the possible excited 
states to density fluctuations and thus ensures the stability of superfluid motion. How- 
ever, as Nozieres says in a very lucid review article [41], "that the T = superfluid 
equations merge with those of ordinary hydrodynamics does not alter the fact that real 
understanding must be based on a microscopic description based on long-range order 
induced by a condensate" . 



(75) 



[lu 2 - (el + 2g'n c e p )] cx [u 2 - (e 2 p + 2gn c e p )][l - 2g X °] - 4g 2 2n c e pX °. 



(76) 



Xnn(l,l')=£ / d2 / d3A a (l,2)G a/ ,(2,3)A /J (3,l')+Xnn(l,l / ) 
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8. Hydrodynamic oscillations in a trapped Bose gas 

In the previous sections, we have been discussing collective oscillations of the con- 
densate and non-condensate in the mean-field or collisionless region. We now turn to 
the collision-dominated hydrodynamic region I <C A, where A is the wavelength of the 
collective mode and I is the collisonal mean- free-path of the elementary excitations [14]. 
Experiments in trapped Bose gases are just now starting to probe the hydrodynamic 
two-fluid region but this region is very rich in new physics and deserves careful study. 

It is unfortunate that the collisionless region at T — described by the Gross- 
Pitaevskii equation is also commonly referred to as "hydrodynamic" theory [1, 2]. I 
guess we are stuck with this terminology, although the terms "quantum hydrodynamics" 
or "superfluid hydrodynamics" would be useful compromises. 

We will limit our analysis to finite temperatures T ~ Tbec, where the number of 
atoms in the non-condensate is comparable to the condensate (N ~ N c ). We will con- 
centrate on the derivation of two-fluid hydrodynamic equations starting from a simple 
microscopic model and try to clarify the physics. These two-fluid hydrodynamic equa- 
tions are expressed in terms of fluctuations of the two components: n c (r, t), v s (r, t) and 
n(r, t), v„(r, t). We can use them to "derive" the two-fluid equations in the conventional 
Landau form, which are written in terms of the fluctuations of local thermodynamic vari- 
ables (temperature, entropy, pressure, etc) . We mainly consider the linearized version of 
the ZGN-type two-fluid equations, which give the hydrodynamic normal modes of the 
coupled system of two components. As a concrete example, we discuss the in-phase rigid 
motion of the equilibrium density profiles of the condensate and non-condensate. 

We need to find equations of motion for both the condensate and the non-condensate. 
As a first approximation [7] , we can use the time-dependent HF Popov equation of motion 
(see Section 6), for the order parameter $(r, t) : 



(77) 



.<9$(r, t) 



dt 



2 



~ + U cxt (r) + gn c (r, t) + 2gn(r, t) 
2m 



$(r,i), 



where n c (r,t) = |$(r, t)\ 2 and h(r,t) is the local non-condensate density. It is clear that 
this equation of motion for n c (r, t) is not closed since it involves 

(78) fi(r,t)=fio(r)+Jn(r,t). 

If we treat the non-condensate in (77) statically, this means that the condensate moves 
in the static Hartree-Fock mean field of non-condensate given by 2gh a (r) [21] More 
generally, we need to have an equation of motion for the fluctuations in Sh(r, t) of the 
non-condensate (ie, the density of the excited atoms), which will be discussed in Section 9. 

We can rewrite equation of motion in (77) in the quantum hydrodynamic variables for 
phase and amplitude, defined by 



(79) 
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where the superfluid velocity is given by 
(80) 



m 



Using (79), (77) can be shown to be equivalent to (all quantities depend on r and t) 



(81) 



m 



+ V • (n c v s ) = 

dv s 



dt 



+ (v s • V)v s 



= -V/i c 



where the superfluid chemical potential is given by (in the dynamic HFP approximation) 



(82) 



Mc(r,t) 



2m, /n 



+ f7 ext (r) + sn c (r, i) + 2gn(r, i). 



The first term in (82) is the "quantum pressure" term, which is ignored in the so-called 
Thomas- Fermi (TF) approximation to the condensate equation of motion [1, 2]. If we 
omit the non-condensate term, equations (81) and (82) describe the well-known Gross- 
Pitaevskii approximation [1]. 

9. Dynamics of the non-condensate atoms 

Following Ref . [7] , we base our discussion on the simplest possible kinetic equation for 
the distribution function /(r, p, t) valid in the semi-classical limit: 



(83) 



IzbT Tiujq (ujq — trap frequency) 
k B T > gn. 



This distribution function is determined by the quantum Boltzmann equation 

df(v,p,t) 



+ Z . Vr - VU(r,t)- V p 
at m 



/(r,P,t) 



dt 



collision 



(84) 

where the effective potential 

(85) U(t, t) = C/ cx t(r) + 2 5 [n c (r, t) + n(r, t)}. 



includes the dynamic HF field. Once we solve (84) for /(r, p, t), we can find the non- 
condensate density from 



(86) 



" (r '* )= /(2^ /(r ' P '* ) - 
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In our simple model, the excited atoms are still particle-like, ie, 
(87) E p (r,t) = ^ + U(v,t). 

This last fact leads to many simplifications but restricts the results to finite temperatures 
close to Tbec- 

We now give the explicit form for the collision term in (84): 



= C 22 [f]+C 12 [f}. 

collision 



The term C 22 [/] describes collisions between excited atoms and is given by [42] , 

n m 1 2 f dpi d p 2 dp 3 3 

CM=2g J — — — (2.) 5(p + Pl -p 2 -p3) 

(89) x 2ttS(E + E 1 -E 2 - E 3 ) [(1 + /)(1 + h)f 2 f 3 - + / 2 )(1 + /s)] • 

In contrast, Ci2[/] describes collisions involving one condensate atom. This kind of 
collision brings about equilibrium between the excited atoms and condensate atoms and 
it is given by (these collision terms were ignored in Ref. [7] ) 

n m o 2 f d Pi dp 2 dp 3 3 

x 2tt5(£' + e c - E 2 - E 3 )(2Trf[S{p - Pl ) - 5(p - p 2 ) - S(p - p 3 )] 

(90) x [n c (l + /!)/ 2 / 3 - n c /i(l + / 2 )(1 + / 3 )] , 

2 

with Ei(r,t) = ^ + U(r,t) and /j = f(r,p i ,t). These two collision terms were first 
derived by Kirkpatrick and Dorfman [43], who considered a uniform gas and worked in 
a frame in which v s (r, t) = 0. Our present calculation [44] is for trapped Bose gas and 
works in the lab frame. As a result, (90) takes into account that a condensate atom has 
energy e c = /j, c + \ravl and momentum p s = mv s , where [i c is defined in (82). We note 
that: 

• C 2 2[/] conserves number, momentum and energy of the colliding excited atoms and 
therefore 

JdpC 22 [f}=0 
JdppC 22 [f} = 

(91) J dpE p C 22 [f] = 

• Ci 2 [/] conserves momentum and energy of the colliding atoms, with the condensate 
atoms having energy e c and momentum rav s . Therefore, we have 
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(92) 



/ dp(p-mv s )C 12 [/] =0 
/ dp(E p - e c )C 12 [f] = 0, 



but 



(93) 



J d P c 12 [/] = r 12 [/]^o. 



The fact that Ti2 is finite follows since C\ 2 does not conserve the number of excited 
atoms; Ti 2 will be referred to as a source term. 

As usual in dealing with the collision-dominated hydrodynamic region, we will assume 
that the collisions C 22 [f] produce local thermal equilibrium among the excited atoms. 
This is described by the local Bose distribution: 



(94) 



/o(r,p,i) 



1 



+t/(r,t)-£] 



where now (3,v n ,U and ft all depend on r, t. One can easily checks that C22[/o] vanishes 
exactly. This uses the key identity for Bose distribution given by (46). In contrast, one 
finds that substituting / = /o in (90) gives 

^ r,, „ 2 f d Pl d P2 dp 3 3 

Ci2[fo}=2 9 J (2n) <5(p 1+ p s -p 2 -p 3 ) 



(95) 



(2tt) 3 (2tt) 3 (2^) a 
x 2w5(E 1 + e c -E 2 - E 3 )(2tt) 3 [S( P - Pl ) - 5(p - p 2 ) - <5(p - p 3 ) 

-m'-»c)] (i + / 01 )/ 02 / 03 , 



x n, 



n c e 



where fi' = fi — ^m(v„ — v s ) 2 ~ fi in a linearized theory. We note that Ci 2 [/o] in (95) 
vanishes when 



(96) 



ie, only if the condensate atoms are in diffusive equilibrium with the non-condensate 
atoms. The same kind of factor arises in the theory of the growth dynamics of a conden- 
sate in trapped Bose gases developed by Gardiner, Zoller and coworkers [45, 46]. 

We next derive hydrodynamic equations for the non-condensate atoms by taking mo- 
ments of the kinetic equation: 



(97) 



- + £ • V r - Vtf(r.t)- V p 
at m 



f(r,p,t)=C 22 [f]+C 12 [f}. 



We obtain [7, 44] 

J dp [kinetic eq.] 



dn 
~dt 



V • nv„ = Tj 
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J (fpp[kinetic eq.] — ► mn ^J^- + v„ • v„ = — VP — nVU — m(v„ — v s )r 

/<9P ~ 2 - 

dpp 2 [kinetic eq.] -» — + V • (Pv„) = — -PV • v n 

(98) +3^c + -m(v„-v s ) 2 -C/]r 12 . 

Here Ti 2 is denned in (93) and 

" (r '* } = / (2^3 /0(r ' P,t) = X3' 93 / 2(Z) 

(") ^ = / (^3^ Mr,p,t)\ Vn=0 = ±±g 6/ 2(z), 

with the local fugacity z(r, t) — e^^ u \ 



10. The two-fluid hydrodynamic equations 

In Sections 8 and 9, we derived a set of 7 equations for non-condensate and the 4 
equations for the condensate. This gives a coupled set of equations for the variables: 

(100) n c ,n, v„,v s ,P and T 12 = J C 12 [fo]- 

^From (95), one sees T12 is function of fo, which in turn depends on /3,v n ,n,n c and jl. 
Thus our system of hydrodynamic equations is closed and can be solved since it involves 
11 local variables and 11 equations. 

For reasons discussed earlier, the explicit effect of the C22 collisions has disappeared 
in the hydrodynamic equations (98). However, the C22 collisions justify and enforce the 
local equilibrium form (94) used for /(r, p,t), which describes the excited atoms. In 
contrast, C12 is still explicitly present in (98) through the Ti2 source term. This means 
that the condensate is not in diffusive equilibrium with non-condensate. We recall from 
(95) that 



(101) r 12 oc 5 2 n c [l-e-' 3 ^)], 

and thus ri2[/o] vanishes only when jl = /i c . One might expect that this condition would 
arise in the case of strong C12 collisions. A detailed study [44] shows that even in this 
case, one must allow for fluctuations in the variable fidiff = A — and thus we cannot 
simply set Ti 2 = 0. Since Ti 2 ^ in (98), we have to generalize the HF Popov description 
to ensure that [compare with the result in (81)] 

(102) ^£ + V -(n c v s ) = -r 12 . 
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This is needed for the continuity equation to be correct 

(103) | + V.j = 0, 
where 

n = n + n c 

(104) j = hv n +n c v s . 

As discussed in detail in Ref.[44], a simple microscopic model based on the neglect of 
anomalous pair correlation functions (i/>(r)i/j(r)) leads to both (102) and (98). 

We can linearize our coupled two- fluid hydrodynamic equations given by (81) and (98) 
around the static equilibrium solution, where /i c o = uq. The fluctuations 6n, 5n Cl <5v„, 5v s 
are precisely the quantities which can be directly measured in trapped gases. The only 
static equilibrium functions we need are ho,n c0 and P , which are given by (99), 

(105) Po = 7^3.95/2(^0), n = ^.93/2(20), 



where z — e~^ 9nc °^ and n c0 is given by the static HFP equation (52) for $(r). We note 
that these static thermodynamic properties for a Bose-condensed gas are quite different 
from the expressions first given by Lee and Yang (1958) for a uniform gas [47]. They 
calculated only the first-order interaction corrections to the thermodynamic properties 
of an ideal Bose gas, using 

^° - ^-T3 #5/2 ( Z = 1) - g^oflcr 

Po A o 

(106) h ~ ^3-53/2(2 = 1) = h cr . 

Here h cr is the critical density at the temperature T for an ideal uniform Bose gas. The 
non-perturbative approximation we use for thermodynamic properties [as illustrated by 
(105)] is consistent with the Lee- Yang results when we expand to first order in g [48]. 

We can use our generalized ZGN equations to derive two-fluid equations in the tradi- 
tional form first given by Landau (1941). In linearized form, these are [6, 49] 

dSn „ 

-^ + v.*j = o 

m ^7 = -V<5P-<5nVC/ext(r) 
at 

m— — = -Vd^ c 
at 

(107) ^ = _ v . ( s< y Vn ), 



THEORY OF EXCITATIONS AT FINITE TEMPERATURES 



25 



where we have included an external potential, s(r, t) is the local entropy density and 



Of course, the Landau equations describe the case when the supcrfluid and normal fluid 
are in dynamic local equilibrium with each other. The Landau form of two-fluid equations 
is very natural for superfluid He, for which they were developed. They involve the usual 
thermodynamic variables: pressure, entropy and temperature fluctuations. However, in 
a trapped gas, the more natural variables are those used in a ZGN-type formulation 
[7, 44], namely h,n c ,v n ,v s and Pdiff- 

In the generalized ZGN equations discussed in Ref. [44] , the fluctuations around static 
equilibrium involve fluctuations in t)Y\i, namely [see (101)] 

(109) ST 12 x g 2 n c0 Sp di f f . 

This leads to an additional equation of motion for Spdi f / , with a characteristic relaxation 
time Tp. This gives rise to a new relaxational mode of the condensate. The ZGN limit [7] 
corresponds to ujt^ ^> 1 while the Landau two-fluid limit corresponds to lut^ <C 1. While 
the relaxational mode is not included, the ZGN equations appear to give a reasonable 
first approximation for the other hydrodynamic modes, as shown in the explicit results 
given in Ref. [48]. We note that Landau two-fluid equations given by (107) have been 
recently used by Shenoy and Ho [50] to work out the hydrodynamic normal modes of a 
trapped superfluid gas. 

The Landau two-fluid equations, which assume total local equilibrium, are valid at 
all temperatures. Our present analysis is restricted to finite T ~ Tbec, but could be 
generalized to the case of very low temperatures. At T = 0, of course, the distinction 
between the collisionless and hydrodynamic region is not a useful one (see, moreover, the 
comment at the end of Section 7). However, at T = 0, we have p n = and p s = p and 
then the Landau two-fluid equations formally reduce to hydrodynamic equations for one 
component: 



These equations are the basis of recent work by Pitaevskii and Stringari [51]. In conjuc- 
tion with the local density approximation [1] 



(108) 



Sj = n c0 Sv s + h Q 5v. 
Sn = Sn c + Sh. 



n 



(110) 




(111) 



fj, c (r,t) = no[n(r,t)] + E/ ext (r) 



where po is the chemical potential of a uniform interacting Bose-condensed gas at T = 0, 
such as given by (62), one may use (110) to find collective mode frequencies past the 
MFA (or GP approxiamation) . 
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11. — In-phase oscillation of the condensate and non-condensate 

Two interesting collective modes exhibited by the ZGN equations involve rigid centre- 
of-mass oscillations of the equilibrium density profiles: 

h(r,t) = n Q (r-rt n (t)), rj n (t) = v n (t) 

(112) n c (v,t) =n c0 (r-77 c (t)), r, c (t) = v a (t). 

One finds [7] an in-phase solution (Kohn mode) with v„ = v s , with a frequency given 
by the trap frequency uj a ; and an out-of-phasc solution, with v„ and v s in opposite 
directions. The latter is the analogue of the second sound mode in liquid 4 He. 

First let us work out the static equilibrium solution / e(? (r, p) of (97), which must 
satisfy 

(113) £ • V r / eg (r,p) - VUo(r) ■ V p / e? (r,p) = C[f eq ] = 0. 
The solution is easily verified to be 

(114) /o(r ' p) = TT^rr, - ; — ' 

where (within the Thomas- Fermi approximation for ^ c0 ) 

= Mco = U ext (r) + 2gn (r) + gn c0 (r) 

(115) U Q (r) - U cxt (r) + 2g[n {v) + n c0 (v)]. 

Using (114), we find ho(r) is given by (105), with 

/ o \ 1/2 

(116) A = , - ; z = e Mfi -°- Ua) = e -A>fl»eo . 

\mk B To ) 

We next consider the in-phase solution described by (112), with v s = v„ = 17 (not 
dependent on the position r). This is a solution of ZGN equations with r] a (t) ~ e ^ at , 
where w a is the parabolic trap frequency in a th direction (x,y,z). The proof proceeds 
as follows (in the linearized theory). Given (112), we have 

Sn c (r,t) = -rj(t) ■ Vn 0c (r) 

(117) 5h(r,t) = -r ] (t)-Vho(v), 



and hence 



m — — — — = — V(5^ c = — V[2gSh + gSn c 



dt 

(118) = V (rj ■ [2gVn + gVn c0 }) . 
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This result is equivalent to [using (115)] 

(119) m0 = - V (77 • VC/ cxt ) , 
with 

(120) U cxt = \m{u? x x 2 + co 2 y y 2 + lo 2 z z 2 ). 
The solution of (119) is clearly 

(121) iia = -U%Va, VaW-e™" 1 . 

In a similar way, one can show that the hydrodynamic equations for the non-condensate 
lead to the same equation of motion (119). The rigid motion of the non-condensate den- 
sity profile corresponds to a distribution function given by 

(122) /(r,p,i) = / (r-T7,p-mr7). 

This can be verified to solve the kinetic equation (97), using the fact that C[fo) — 0, and 
leads to the expression in (112), namely 

(123) n(r, t) = j f(r, p, t) = n (r - V ). 

This in- phase dipole mode was first exhibited by ZGN in the hydrodynamic region [7]. 
It is an important test of the proper treatment of the non-condensate dynamics. This 
kind of mode is generic, existing in the collisionless region as well [29, 34]. 

12. First and second sound in a uniform Bose-condensed gas 

The simplicity of the ZGN two-fluid equations allow a very transparent discussion of 
first and second sound in a uniform gas [48]. We find two sound-wave solutions: to = uiq, 
to = u 2 q. First sound has a velocity given by (to lowest order in g) 

(124) „? *4_ + ?£no 

3 mno m 

and mainly involves the non-condensate, with v n 3> v s (but in- phase). Second sound has 
a velocity given by (to lowest order in g) 

(125) ul = ^o, 

and mainly involves the oscillation of the condensate, with v s 3> v n (and out-of-phase). 
It is clear that second sound is the collision-dominated hydrodynamic mode which is the 
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analogue of the collisionless Bogoliubov phonon (42) which wc discussed in Section 4. The 
frequency vanishes (becomes "soft" ) above Tbec and it is the expected Goldstone mode 
of the Bose broken-symmetry [16, 17]. In contrast to liquid 4 He, both second and first 
sound in a Bose gas involve density fluctuations (ie, second sound is not a temperature 
wave in a gas!). 

Can we study first and second sound pulse propagation in cigar-shaped traps? In 
this volume, Kctterlee [31] reviews the beautiful pulse propagation experiments carried 
out at MIT at low temperatures in the collisionless region [52] . If we could get into the 
hydrodynamic region, one should be able to see two pulses propagate — roughly with the 
speeds u\ and u-i given above. This would be very dramatic and a direct confirmation of 
the superfluid dynamics described by the two-fluid equations. Inside the broad, relatively 
uniform condensate distribution along the axis of the cigar-shaped trap, one could excite 
both first and second sound pulses, with the relative intensity of each propagating pulse 
given by [53] 

(126) ^5{z- Ui t). 

The velocity of second sound would give [using (125)] a direct measurement of the con- 
densate density n c0 . Perhaps of more interest, the velocity of the first sound pulse would 
give a relatively direct measurement [using (124)] of the non-condensate density n n "un- 
derneath" the condensate. 

13. Concluding remarks 

In Sections 9 and 10, we discussed the linearized form of the generalized two-fluid 
equations [44] which describe how the condensate relaxes to equilibrium with the excited 
atoms. These hydrodynamic equations with T\2 ^ (i.e. jl ^ fi c ) should also be useful in 
the study of how a condensate grows and comes into equilibrium with the thermal cloud 
[31]. While based on a different formalism, our theory has points of contact with the 
work by Gardiner, Zoller and co-workers on the kinetics of condensate growth [45, 46]. 
However, our analysis, while assuming the existence of a condensate, takes the dynamical 
evolution of the non-condensate fully into account. The two- fluid dynamics of trapped 
Bose-condensed gases promises to be a very rich subject, as discussed in Ref. [44] . The 
region where Ti2 ^ is not easily studied in superfluid 4 He, largely because local thermal 
equilibrium is almost too easy to reach in a liquid. However, some years ago Pitaevskii 
[54] treated the regime very close to T\, when the superfluid density p s (r,t) was not 
necessarily in equilibrium with the normal fluid density Pn(r, t). 

One of the interesting aspects of trapped Bose gases is that even if the centre of the 
trap is at large enough density so that one is in the hydrodynamic region described by 
dynamic local thermal equilibrium, this description will always break down in the low 
density non-condensate thermal cloud far enough from the trap centre. Recent theoretical 
studies [55, 56] have shown that the cross-over region between the hydrodynamic and 
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collisionless regions in this low density region plays the dominant role as the source of 
damping of certain kinds of hydrodynamic oscillations. 

More generally I would like to emphasize that the study of collective modes in trapped 
Bose gases finally provides us with the opportunity of making a quantitative test of the 
complex dynamics of a system with a Bose-broken symmetry. As I discuss elsewhere in 
this volume [9], the complexity of dealing with a liquid such as superfluid 4 He never al- 
lowed one to make a thorough confrontation of the microscopic theory of Bose-condcnscd 
systems with experiment [11]. Bose gases finally allow us to do this, in a much cleaner 
fashion [31, 57]. I would hope that future studies on trapped Bose gases will also stimulate 
new interest in the dynamics of superfluid 4 He as a Bose-condensed liquid. 

* * * 

The work presented here on the two-fluid hydrodynamic equations has been done 
in close collaboration with Eugene Zaremba and Tetsuro Nikuni. In addition, I have 
benefitted from stimulating discussions with many of my colleagues, and especially would 
like to thank K.Burnctt, M.L. Chiofalo, E. Cornell, W. Ketterlc, L.P. Pitacvskii, G.V. 
Shlyapnikov, D.M. Stampcr-Kurn, S.Stringari and P. Zoller. 

These lectures are partly based on research done while participating in the 1998 BEC 
Workshop sponsored by the Institute of Theoretical Physics in Santa Barbara. 

My research is supported by NSERC of Canada. 
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Figure Captions 

1. Hartree-Fock self-energies in a normal Bose gas. 

2. Various interaction terms involving condensate (wiggly line) and non-condensate 
(solid line) atoms. 

3. Self-energy diagrams in the Bogoliubov approximation (T = 0). 

4. Self-energy diagrams in the first-order Popov approximation. The propagators are 
for an ideal Bose gas. 

5. Self-consistent Hartree-Fock-Bogoliubov (HFB) approximation for self-energies. 

6. Ladder diagram approximation for the i-matrix. 



